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Emission and propagation of Saturn kilometric radiation: 
magneto-ionic modes, beaming pattern and polarization state 

L. Lamy 1,2 , B. Cecconi 1 , P. Zarka 1 , P. Canu 3 , P. Schippers 4 , W. S. Kurth 4 , R. L. 
Mutel 4 , D. A. Gurnett 4 , D. Menietti 4 , P. Louarn 5 

Abstract. The Cassini mission crossed the source region of the Saturn kilometric ra- 
diation (SKR) on 17 October 2008. On this occasion, the Radio and Plasma Wave Sci- 
ence (RPWS) experiment detected both local and distant radio sources, while plasma 
parameters were measured in situ by the magnetometer (MAG) and the Cassini Plasma 
Spectrometer (CAPS). A goniopolarimetric inversion was applied to RPWS 3-antenna 
electric measurements to determine the wave vector k and the complete state of polar- 
ization of detected waves. We identify broadband extraordinary (X) as well as narrow- 
band ordinary (O) mode SKR at low frequencies. Within the source region, SKR is emit- 
ted just above the X mode cutoff frequency in a hot plasma, with a typical electron-to- 
wave energy conversion efficiency of <~ 1% (2% peak). The knowledge of the k- vector 
is then used to derive the locus of SKR sources in the kronian magnetosphere, that shows 
X and O components emanating from the same regions. We also compute the associated 
beaming angle at the source 9' = (k,-B) either from (i) in situ measurements or a model 
of the magnetic field vector (for local to distant sources) or from (ii) polarization mea- 
surements (for local sources). Obtained results, similar for both modes, suggest quasi- 
perpendicular emission for local sources, whereas the beaming pattern of distant sources 
appears as a hollow cone with a frequency-dependent constant aperture angle: 8' = 75°± 
15° below 300 kHz, decreasing at higher frequencies to reach 6»'(1000 kHz) = 50°±25°. 
Finally, we investigate quantitatively the SKR polarization state, observed to be strongly 
elliptical at the source, and quasi-purely circular for sources located beyond approximately 
2 kronian radii. We show that conditions of weak mode coupling are achieved along the 
ray path, under which the magneto-ionic theory satisfactorily describes the evolution of 
the observed polarization. These results are analyzed comparatively with the Auroral Kilo- 
metric Radiation (AKR) at Earth. 



1. Introduction 

Forty years of remote and in situ measurements of 
Earth-orbiting spacecraft led the auroral kilometric radia- 
tion (AKR) to become understood enough to be used as a 
basis for investigating other planetary auroral radio emis- 
sions. AKR is a powerful non-thermal emission with com- 
plex properties generated above the terrestrial atmosphere 
[Benediktov et al, 1965; Gurnett, 1974] by the Cyclotron 
Maser Instability (CMI) [Wu and Lee, 1979; Treumann, 
2006] from unstable hot electrons accelerated at energies of 
a few keV [Louarn et al., 1990; Ergun et al, 2000]. 

Saturn kilometric radiation (SKR) is the kronian equiva- 
lent of AKR, knowledge of which essentially lies on remote 
observations of Voyager (flybys in 1980 and 1981) [Kaiser 
et al, 1984], Ulysses (distant observations in the 1990s) and 
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Cassini (in orbit since mid-2004) [Kurth et al., 2009]. Its 
spectrum extends from a few kHz to 1200 kHz, for altitudes 
between 5 and 0.1 Rs (1 Rs = 1 Saturn radius = 60268 km). 
Like AKR, SKR displays remote properties consistent with 
the CMI mechanism [Zarka, 1998]. It is emitted mainly in 
the X mode, with fainter O mode [Lamy et al., 2008a; Cec- 
coni et al, 2009]. Radio sources lie along magnetic field 
lines globally associated with atmospheric aurorae [Lamy 
et al, 2009], with correlated power radiated in radio and 
UV ranges, and a conversion efficiency between precipitated 
energy and radiated SKR around 0.5% [Kurth et al, 2005; 
Clarke et al., 2009]. 

The apparent beaming pattern of SKR has been esti- 
mated with two different approaches. On the one hand, di- 
rect measurements showed a frequency-dependent beaming 
angle, varying within the range 40° — 50° (40° — 80° respec- 
tively) in the southern (northern respectively) hemisphere, 
and slightly decreasing with frequency [Cecconi et al, 2009]. 
On the other hand, modeling work of the SKR visibility ob- 
tained best results with a frequency-dependent beaming an- 
gle, varying from ~ 70° at 3 kHz to ~ 60° at 1000 kHz for 
both hemispheres, computed from a CMI-unstable loss cone 
distribution with 20 keV electrons [Lamy et al, 2008b]. 

SKR was first observed to be almost purely circular polar- 
ized, from instantaneous and statistical equatorial and low 
latitudes observations [Ortega-Molina and Lecacheux, 1990; 
Cecconi et al., 2006; Lamy et al, 2008a], but high latitude 
measurements then revealed that kilometric waves become 
strongly clliptically polarized when observed from latitudes 
between 30° and 60° [Fischer et al, 2009]. 

On 17 October (day 291) 2008, the Cassini spacecraft 
crossed the SKR source region for the first identified time 
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at an unusual location near midnight Local Time (LT) at a 
5 Rs distance and latitudes between —60° and —70°, bring- 
ing crucial in situ measurements required for a better un- 
derstanding of the emission [Kurth et al., 2010]. 

On this occasion, we investigated quantitatively the lo- 
cal wave and plasma properties [Lamy et al, 2010a] (here- 
after paper I), briefly summarized below. An overall source 
region, characterized by intense signal close to f ce was de- 
tected by the Radio and Plasma Wave experiment (RPWS) 
between 0812 and 0912 UT. Within this interval, we identi- 
fied 3 unambiguous events (dividing in 8 individual RPWS 
measurements) strictly below f ce , that, according to exten- 
sive AKR studies at Earth, were considered as actual tra- 
versed sources. The footprints of magnetic field lines sup- 
porting local and distant SKR sources revealed a spiral au- 
roral oval starting from nightside very high latitudes, in- 
dicating an enhanced auroral activity, proposed by Bunce 
et al. [2010] to result from a solar wind compression. Like 
at Earth, the auroral plasma is hot and tenuous, but unlike 
Earth, the SKR source region was crossed at much higher 
altitudes (4.1 Rs), where no auroral cavity (at the edges 
of which AKR propagation is affected) was detected. SKR 
frequencies observed below f ce were found compatible with 
CMI resonance frequencies computed with the observed 6 to 
9 keV hot electrons as resonant particles and for amplifica- 
tion perpendicular to the magnetic field. Interestingly, these 
hot electrons display ring or shell-like distributions similar 
to those responsible for AKR generation at Earth [Schippers 
et al, 2010], and have been used to compute CMI growth 
rates high enough to account for observed SKR intensities 
[Mutel et al, 2010]. 

This event is of particular interest, for remote and lo- 
cal SKR sources were measured simultaneously. In this 
study, we extend the analysis of paper I to investigate 
polarization and beaming properties of SKR waves, along 
their propagation from their source region. After having 
described Cassini's instrumentation and observational pa- 
rameters (section 2), as well as useful wave properties in a 
magnetized plasma (section 3), we first identify and char- 
acterize the different magneto-ionic modes of SKR (section 
4). Then, we compare the location of regions from which 
observed modes emanate (section 5) and investigate their 
beaming pattern, derived for local to distant sources (sec- 
tion 6). Finally, we quantify the SKR polarization state 
at the source and its evolution with propagation (section 
7) and relate it to predictions of the magneto-ionic theory. 
These results are compared to other auroral planetary radio 
emissions, especially to the terrestrial case. 



wave, namely its full state of polarization, defined by the 
four Stokes parameters S, Q, U and V [Kraus, 1966], and, 
with the additional hypothesis of transverse electromagnetic 
waves (k.E = 0), the k- vector direction, defined by two an- 
gular coordinates in the spacecraft frame. As the latter is 
determined by the direction perpendicular to the polariza- 
tion ellipse, it cannot be unambiguously derived for a purely 
linearly polarized wave. Out of this limit, a direction yields 
two possible senses of k, among which we select the one 
corresponding to the source the closest to the planet. Each 
quasi-instantaneous measurement at a given frequency thus 
characterizes the radio wave radiated by the most intense 
point source detected at that time and frequency. 

Over the 3.5 to 1500 kHz spectral range of interest 
here, HFR measurements were acquired with logarithmically 
spaced frequency channels distributed within three consec- 
utive bands from 3.5 to 325 kHz, with a spectral resolution 
Sf/f = 5%, and linearly spaced channels within the high 
frequency band HF1 above 325 kHz, with a fixed resolution 
of 25 kHz. 

In the supplementary material of paper I, we described 
how radio data were specifically processed and selected along 
this interval to obtain reliable measurements. We remind be- 
low selection criteria applied to the degree of circular polar- 
ization V, the signal-to-noise ratio SNR, and the z r param- 
eter quantifying the variation of the autocorrelation on the 
w monopole between two consecutive 2-antenna measure- 
ments: 0.05 < |V| < 1.1 (to remove unpolarized or aberrant 
measurements), SNR > 20 dB simultaneously on each an- 
tenna of the full 3-antenna measurement, |Vj > 0.7 when 
20 < SNR < 45 dB (to remove specifically weakly circularly 
polarized emissions with high SNR as low frequency vari- 
able background and narrowband emissions), and \z r \ < 0.05 
(corresponding to a tolerance of signal variation of 10 % be- 
tween 2 consecutive sets of measurements) . Out of Figures 1 
and 2 that used all available data excluding radio frequency 
interference, all other Figures were built from the above data 
selection. 

The RPWS experiment also includes three orthogonal 
search coil magnetic antennas aligned with the x, y and 
z axis of the spacecraft, and a Medium Frequency Receiver 
(MFR) [Gurnett et al, 2004], in charge of measuring the 
orthogonal magnetic components of electromagnetic waves. 
Along the present interval, the MFR provided measurements 
of the B z spectral density between 0.02 kHz and 12 kHz, 
with logarithmically spaced frequency channels and a spec- 
tral resolution of 7%. 



Finally, plasma parameters were obtained along the 
2. Instrumentation and observational paramefceESscraft trajectory thanks to simultaneous magnetic field 

and electrons in situ measurements. High resolution ob- 



Together with the three electric antennas u, v and w, 
the RPWS experiment contains a High Frequency Receiver 
(HFR) which measures the wave electric power spectral den- 
sity between 3.5 kHz and 16.125 MHz, including the usual 
spectral range of SKR waves [Gurnett et al, 2004]. The 
HFR records a set of two auto-correlations and a complex 
cross-correlation of input signals sensed on a pair of anten- 
nas, hereafter simply called a 2-antenna measurement. 

During the time interval investigated in this study, the 
HFR was set up in the 3-antenna operating mode, consisting 
of two consecutive quasi-instantaneous 2-antenna measure- 
ments over the (u, w) and (v, w) pairs of monopoles. This 
simulates a real 3-antenna observation when the observed 
signal does not vary significantly between both 2-antenna 
measurements. 

Under the point source assumption, a goniopolarimetric 
(GP) inversion (compounding the words goniometric - for 
angles - and polarimetric - for polarization - analysis) ap- 
plied to each 3-antenna measurement [Cecconi and Zarka, 
2005] derives the six physical parameters of the observed 



servations of the magnetometer (MAG) [Dougherty et al., 
2004] provided the electron cyclotron frequency f ce and the 
local magnetic field vector B, while electron moments were 
computed from measurements of the Electron Spectrometer 
(ELS) of the Cassini Plasma Spectrometer (CAPS) [Young 
et al., 2004], giving the plasma frequency f pe as well as the 
temperature and the total energy carried by hot electrons. 

3. Waves in a magnetized plasma 

In this section, we recall some useful wave properties in 
a tenuous and magnetized plasma (f pe <g f ce ), as the one 
observed in the kronian auroral region. 

3.1. Cold plasma dispersion and magneto-ionic 
theory 

It is well known that electromagnetic waves can prop- 
agate in a magnetized plasma in several characteristic (or 
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natural) modes, defined by indices of refraction N which 
are solutions of the dispersion equation. For a homogeneous 
cold plasma, neglecting the ions motion (high frequency ap- 
proximation) and collisions, the dispersion relation, known 
as the Altar- Appleton-Hartree or Appleton-Lassen equation 
[Appleton, 1932; Lassen, 1927], is: 



N 



1 - 



2X(1-X) 



2(1 -A)- Y 2 sin 2 (9=FA 
with A = \Jy a sin 4 9 + 4Y 2 (1 - X) 2 cos 2 9 



(1) 



where X = (f P e/f) 2 and Y — fee/ f are the characteristic 
frequencies of the medium where / is the wave frequency, 
6 = (k, B) is the wave propagation angle, and the signs — 
and + refer to the extraordinary and ordinary solutions, re- 
spectively. This expression is the basis of the magneto-ionic 
theory. 

In the general case of oblique propagation, the upper 
branch of the extraordinary solution and the ordinary solu- 
tion are often labelled R-X and L-O modes, as an interme- 
diate case between R,L modes (solutions for parallel propa- 
gation) and X,0 ones (solutions for perpendicular propaga- 
tion) [Goertz and Strangeway, 1995]. As the lower branch 
of the extraordinary solution, namely the Z mode, also dis- 
plays R-X characteristics above f pe (and L-X below), we 
will hereafter use the simple X, O and Z denomination for 
clarity. 

Superluminous (N < 1) X and O modes are defined 
above their cutoff frequency, fx ~ /ce(l + (fpe/ fee) 2 ) and 
So — fpe respectively, and tend asymptotically toward a 
light wave at high frequencies (N — 1), so that they can 
ultimately propagate freely through space. In contrast, the 
Z mode, subluminous above f pe and superluminous below, 
is trapped between its cutoff frequency fz ~ fpe/ fee and the 
upper hybrid (resonance) frequency fun = {fpe + fce) 1 ^ 2 ~ 
f cc (l + l/2(/ pe / fee) 2 )- As fuH is always lower than fx, a 
frequency gap separates X and Z components. 

In addition to the refractive index, an important char- 
acteristic of magneto-ionic modes given by the dispersion 
equation is the wave polarization, which is elliptical in the 
general case [Melrose, 1980; Budden, 1985]. The transverse 
component of the polarization, defined as the axial ratio T 
of the polarization ellipse in the plane perpendicular to k, 
is related to the refractive index by: 



Ni 



T x o - Ycos9 



T can be expressed as a function of X, Y and 9: 



T xo = A^sign(A)\/l+A 2 



(2) 



(3) 



with A - 



Y sin 2 9 
2(1 — X) cos 9 



where the signs — and + refer to the extraordinary and or- 
dinary solutions, respectively. The longitudinal component 
of the wave polarization (along k) is neglected, being of the 
order of XY. The solutions of equation 3 are orthogonal 
(T X T C = -1). 

The axial ratio is related to the Stokes parameters [Kraus, 
1966], and the normalized degree of circular polarization V 
is given by: 

9T 

v = YTt- 2 (4) 

We retrieve the limiting cases of purely circularly polar- 
ized waves for parallel propagation (T = ±1, V = ±1), and 
purely linearly polarized ones for perpendicular propagation 



(T = 0, ±oo, V = 0), whereas the wave polarization is ellip- 
tical in the general case. 

The sign of V and T indicates the sense of polarization: 
for A > 0, T x ,V X < corresponds to a RH polarized wave, 
and T Q ,V > to a LH polarized wave, while the opposite 
situation occurs for A < 0. Thus, for 9 > ir/2, which corre- 
sponds to the general case of SKR emission in the southern 
kronian hemisphere, X and subluminous Z mode are LH 
polarized while O and superluminous Z mode are RH polar- 
ized. 

3.2. Hot plasma dispersion 

Similarly to the terrestrial auroral cavities, the kronian 
auroral region is dominated by weakly relativistic (hot) elec- 
trons while thermal (cold) electrons are typically 5 times less 
dense. The cold plasma approximation is consequently very 
rough and relativistic effects need to be taken into account. 
In the frame of AKR generation, several authors solved the 
wave dispersion equation in a hot plasma, and showed that 
this essentially results in lower characteristic frequencies of 
the extraordinary branches, which in turn favors CMI reso- 
nance close to fee [Wu et al, 1981; Wong et al, 1982; Pritch- 
ett, 1984; Le Queau et al., 1984a, b; Winglee, 1985; Louarn 
and Le Queau, 1996a]. 

For the purpose of this article, we are interested in the 
modified X mode cutoff, that was expressed from various 
typical hot electron distributions as: 



/4i + «(^) 2 + /3( 

Jce 



Ve ,2, 



(5) 



where v e is the characteristic velocity of dominant hot 
electrons, c is the speed of light and (a,/3) are coefficients 
related to the chosen distribution function. Winglee [1985] 
computed fx for a Maxwellian (a = l,/3 = —5/2) and a 
Dory, Guest and Harris (DGH) distribution at first order 
(a — I, /3 — —9/2) under the condition (v e /c) 2 < (f pe /fce) 2 . 
Louarn and Le Queau [1996a] derived a more general ex- 
pression of fx for an idealized ring-like distribution (a — 
1/2, p = —1/2), and underlined that more realistic electron 
distributions would not affect it significantly. Equation 5 
shows that fx can be lower than f ce for energetic enough 
electrons, in agreement with AKR observations. 

4. SKR modes of emission and propagation 

RPWS observations of day 291, 2008 are displayed in Fig- 
ure 1, between 0600 and 1100 UT, in the form of three dy- 
namic spectra, bringing complementary informations on the 
nature of Saturn kilometric radiation. Figure la displays 
the spectral power density recorded on the w monopole be- 
tween 3.5 and 1500 kHz. While intense electromagnetic nar- 
rowband (NB) emissions are seen around 5 kHz [Ye et al., 
2010], together with weak sporadic broadband emission be- 
low 10 kHz related to auroral hiss [Kopf et al, 2010], SKR 
is dominant between ~ / ce and 1000 kHz. As identified in 
paper I, the overall source region corresponds to enhanced 
SKR close to f ce between 0812 and 0912 UT (vertical dashed 
lines), with three successive events A, B and C measured at 
frequencies below f ce (delimited by vertical dotted-dashed 
lines). 

4.1. Magneto-ionic modes 

Similarly to AKR, remote observations identified SKR 
propagating in both free-space X and O modes, but unlike 
AKR, no Z mode SKR has been detected so far. 

Here, we can directly identify magneto-ionic modes in 
Figure lb, that displays a dynamic spectrum of circular po- 
larization, thanks to wave characteristics (frequency domain 
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Figure 1. RPWS-HFR dynamic spectra of (a) electric spectral density recorded by the w monopole, (b) 
degree of circular polarization V, as derived from a 3-antenna GP inversion, over day 291 of year 2008 
between 0600 and 1100 UT. (c) RPWS-MFR dynamic spectrum of magnetic spectral density recorded 
simultaneously along the z axis of the spacecraft, between 3 and 12 kHz. A daily background, computed 
for each frequency, has been subtracted to individual time-frequency measurements. White dashed and 
solid curves indicate f ce and f pe , as derived from MAG and CAPS measurements. Vertical dashed 
lines remind the source region identified in paper I, where SKR is enhanced close to / ce , while vertical 
dotted-dashed lines delimit events A, as well as B and C, for which / < f CE . 



and sense of polarization) discussed in section 3. LH polar- 
ized emission (black, V > 0), dominant over the whole SKR 
spectrum and continuously observed from local (/ ~ f ce ) to 
distant (/ > f ce ) sources, is consistent with predominant X 
mode radiated and propagating at 9 > n/2 (approximately 
outward the planet), according to equation 4. At Earth, 
Z mode AKR corresponds to broadband weak emission [de 
Feraudy et al, 1987, and references therein] that can extend 
significantly below f ce , which is not observed here. However, 
the HFR spectral resolution (5f = 450 Hz at 9 kHz) is not 
sufficient to resolve the forbidden frequency gap between X 
and Z modes (fx — fuH ~ 45 Hz, using f pe = 900 Hz as de- 
rived in paper I) and, therefore, we cannot formally exclude 
the presence of the latter at the lower edge of the X mode 
emission. 

SKR also displays sporadic RH emission (white, V < 0) 
below 40 kHz. This corresponds either to O mode emitted 
in the same direction as the main LH component (9 > 7r/2), 
or to X mode radiated in the opposite direction. Indeed, at 



the traversed source, VB lies at a few degrees relative to B, 
allowing an X mode propagation window of a few degrees 
toward the planet. Here, we favor O mode because this 
emission is fainter than the predominant LH X mode, and 
appears at the lower edge of its envelope, as expected from 
fo < fx, and in agreement with previous observations at 
Saturn [Ceccom et al, 2009] or the Earth [Panchenko et al., 
2008]. For comparable frequencies between 15 and 30 kHz, 
the X/O intensity ratio reaches ~ 10 2 . The overall SKR pat- 
tern displays highly circularly polarized emission (|Vj ~ 1), 
except around / ce , where V decreases while the emission 
level remains high. This is a consequence of elliptical po- 
larization observed for local sources, that is investigated in 
details in section 7. 

Finally, we note that low frequency NB emissions and 
auroral hiss display clear Z mode characteristics, with LH 
polarization between f pe (< 2 kHz after 0700 UT) and f ce , 
and RH polarization below f pe , indicating wave propagation 
at 6 > n/2. 
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Figure 2. Characteristic frequencies between 0812 and 0912 UT of day 2008-291, expressed relatively 
to f ce . The SKR cutoff frequency f cut , determined with an accuracy of 0.5 %, is displayed in black, and 
fee, in dashed blue. The curves labelled (1), (2) and (3) plot the X mode cutoff frequency fx computed 
from (1) maxwellian, (2) DGH [Winglee, 1985] and (3) ring-like electron distributions [Louarn and Le 
Queau, 1996a], using fee, f pe and the kinetic energy of (6 to 9 keV) hot electrons, as derived from MAG 
observations and CAPS-ELS 200 s-averaged moments. 



4.2. Electromagnetic waves 

Figure lc shows a dynamic spectrum of the wave mag- 
netic field component measured along the z axis of the space- 
craft between 3.5 and 12 kHz. In spite of a lower sensitivity 
than electric measurements of Figure la, a clear signal is 
detected within the SKR envelope, as well as a fainter one 
within the trace of NB emission at the end of the inter- 
val, giving a direct evidence of the electromagnetic nature 
of these emissions. In addition, Figure lc confirms the ob- 
servation of SKR events / < / ce at timings consistent with 
the ones identified in paper I. 

Due to the uncertainty on the wave direction with respect 
to the z axis and the lack of high SNR measurements, it was 
not possible to derive reliable, accurate values of the refrac- 
tive index. However, the ratio of the wave electric field E 
and c.B was often close to unity at frequencies slightly above 
f ce , consistent with N ~ 1. 

4.3. X mode cutoff 

Based on remarkably similar AKR and SKR features, the 
CMI has been proposed as a common mechanism respon- 
sible for both radiations. Through this process, resonant 
electrons can amplify a background of radio waves (see for 
instance [Zarka et al, 1986]) at frequencies close to / ce , pri- 
marily in the X mode. To check the robustness of this in- 
terpretation, we want to verify that the SKR local emission 
frequency remains higher than the theoretical X mode cutoff 
frequency fx- 

The SKR low frequency cutoff f cu t , as determined by the 
lower limit of the SKR envelope, is plotted in black in Fig- 
ure 2 over the source region and relatively to f ce . In paper 
I, we showed that the CMI resonance frequency matches 
fcut only for crossed sources A, B and C, distributed in 8 
measurements f cu t < fee- 

From in situ measurements of / ce , f pe and v e , we calcu- 
lated the three types of fx given in section 3.2 along the 
Cassini trajectory, shown in Figure 2 by curves (1), (2) 
and (3). X mode cutoffs from Winglee [1985] are always 
lower than the CMI frequency. Nevertheless, as the condi- 
tion (v e /c) 2 < (fpe/ fee) 2 is poorly satisfied in the source 
region, the real X mode cutoff should lie above curves (1) 
and (2). A more realistic cutoff, derived by Louarn and Le 
Queau [1996a], is given by curve (3). It matches almost all 
feut of events A, B and C within error bars, with slightly 
higher average values. However, no uncertainty was consid- 
ered on fx since it mainly depends on the technique used 
to derive the temperature of hot electrons. Whereas the 
electron moments used here were computed through a di- 
rect integration of CAPS-ELS data, averaged over 200 s to 
obtain a continuous series, a forward modeling of individ- 
ual CAPS-ELS spectra [Schippers et al, 2008] leads to less 



continuous but somewhat hotter temperatures (up to 1 keV 
more), which in turn result in ~0.02% lower fx, in better 
agreement with SKR cutoff. 

Investigating the regimes authorized by the CMI reso- 
nance, Louarn and Le Queau [1996a] derived a critical elec- 
tron kinetic energy E cr = (3/8)(/ pe // ce ) 2 m c c 2 , above which 
X mode is favored over Z mode, and conversely below. In pa- 
per I, we computed ratios f ve / fee = 0.05 — 0.09 within the 
source region, corresponding to E cr = 0.5 — 1.5 keV, well 
below the 6 — 9 keV range of dominant (and resonant) hot 
electrons. Finally, O mode growth rates have been shown to 
be significantly lower than the growth rate of dominant (X 
or Z) mode, whatever the f pe /fce ratio [Wu and Lee, 1979; 
Winglee, 1985]. These results are fully consistent with the 
SKR modes identified above. 

4.4. Electron-to-wave energy conversion efficiency 

The energy conversion efficiency from unstable electrons 
to amplified radio waves has been estimated in previous 
studies by two methods. 

A macroscopic approach consists of comparing the to- 
tal radio power radiated over its full spectrum assuming a 
known solid angle, to the total input (or dissipated) electron 
energy, assuming that the latter is a good approximation 
of the free energy available to the CMI resonance. Doing 
so, previous estimations led to typical efficiencies of 1% for 
AKR [Gurnett, 1974] and 0.5% for SKR [Kurth et al, 2005]. 
However, the input energy only includes precipitating elec- 
trons, while CMI-unstable electron distributions found in 
radio source regions involve mirrored particles, with maxi- 
mum growth rates for electrons at large pitch angles, namely 
close to their mirror point. Thus, the input energy only gives 
a rough estimate of the available free energy. 

Alternately, a microscopic approach directly compares ra- 
dio and particle measurements within the source region us- 
ing a few assumptions. At Earth, Benson and Calvert [1979] 
computed the total AKR power radiated over its full spec- 
trum from the electric field strength measured locally and 
compared it to the local particle flux to derive efficiencies 
between 0.1 and 1%. 

A similar computation can be performed here. In paper 
I, we derived a typical electric field strength of ~ 1 mV.m , 
estimated from the mean spectral power density measured 
on the single w monopole and integrated over a 5 kHz wide 
emission. Here, we are interested in estimating the peak 
Poynting flux S=eocE 2 /2 and electric field amplitude close 
to fee, for the 8 individual measurements of events A, B 
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Figure 3. (a) Location of radio sources as seen from Cassini, plotted with a field of view (FOV) of 
126°, and (b) magnetically projected down to the planet over the time-frequency interval [0812,0912] UT 
(source region) and [/ ce -l kHz,/ cc +40 kHz] (highest resolution). This Figure is identical to Figure 2a, b 
of paper I except that black and blue symbols here distinguish X (RH polarized) from O (LH polar- 
ized) mode emissions. Local X mode sources, corresponding to frequencies between f ce — 1 kHz and 
fee + 1 kHz, are displayed in red. Both X and O mode SKR sources emanate from the same magnetic 
field lines, clustered on the dawn side, and whose footprints form a spiral shape (see also the radio map 
displayed by Figure 2c in paper I). 



Table 1. Table of radio and electron observables observed during events A, B and C. 



SKR sources 


A 


B 


C 


DOY 2008-291 (h:min:s) 


08:38:53 


08:39:25 


08:39:57 


08:40:29 


08:52:45 


08:53:17 


08:55:25 


08:55:57 


Peak Poynting flux (xl0 _y W.m~' z ) 


9.2 


3.5 


9.3 


6.0 


2.3 


1.2 


0.96 


3.6 


Electric field amplitude (mV.m - 1 ) 


2.6 


1.6 


2.6 


2.1 


1.3 


0.96 


0.85 


1.6 


Total radiated power (xl0~ 4 crgs.s~ 1 .cm~ 2 .sr _1 ) 


5.9 


2.0 


5.7 


3.0 


1.2 


0.91 


0.71 


2.7 


E e (> 1 keV) (xlO-- 1 ergs.s-^cm-^.sr" 1 ) 


4.0 


4.0 


3.0 


3.0 


2.2 


2.2 


3.6 


2.1 


Electron-to-wave efficiency (%) 


1.5 


0.50 


1.9 


1.0 


0.55 


0.41 


0.20 


1.3 



and C, once corrected from the geometric projection of the 
wave electric field onto the antennas, and integrated over the 
instrumental bandwidth of each filter (~ 0.5 kHz). There- 
fore, we use the magnitude of the Poynting vector given 
by the 3-antenna GP inversion assuming N = 1. The 
peak Poynting flux for the first measurement of event A 
reaches 9.2 x 10~ 9 W .m 2 , corresponding to an electric field 
strength of 2.6 mV.m , slightly higher than our previous 
estimation. Considering that the signal observed at 10 kHz 
reached the 1% occurrence level, whose typical spectrum 
can be approximated by a constant signal over a 800 kHz 
bandwidth (see Figure 8 of Lamy et al. [2008a]), we com- 
pute the equivalent power flux radiated along an entire field 
line 1.5 x 10 -5 W.m -2 . Then, anticipating results of section 
6, we consider that the SKR is beamed within an emission 
cone of solid angle ~ 0.5 sr. The corrected power flux in 
47r sr then reaches 5.9 x 10 -7 W.m _2 .sr -1 , or equivalently 
5.9 x 10~ 4 ergs.cm _2 .s _1 .sr~ 1 . 

On the other hand, the available free energy can be esti- 
mated by the CAPS-ELS derived energy flux for electrons 
> 1 keV, averaged over all anodes, assuming that hot elec- 
trons are CMI-unstable over the SKR source region. This 
energy flux represents more than 97% of the total electron 
energy flux over the source region, which illustrates that the 
energy flux of cold electrons is negligible. For our example, 
E e (> 1 keV) = 0.040 ergs.s _1 .cm -2 . Compared to the peak 
SKR energy radiated, we obtain a final electron-to-wave en- 
ergy conversion efficiency of 1.5%. 



Table 1 summarizes values obtained for each of the 8 mea- 
surements constituting source events A, B and C. The final 
mean efficiency reaches 0.9% with a standard deviation of 
0.6%, in excellent agreement with previous estimations. 

5. Source location 

Under the assumption of a cyclotron emission at / = f ce 
and straight line propagation, RPWS' ability to retrieve the 
direction of incoming radio waves, coupled to a magnetic 
field model (SPV model [Davis and Smith, 1990] including 
a simple current sheet [Connerney et al., 1983]), enables us 
to derive both the spatial location of radio sources, and the 
beaming angle 6 between the wave vector k and the local 
magnetic field vector B at the source (see [Cecconi et al., 
2009] for details). 

The overall spatial distribution of local and distant SKR 
sources, mixing all polarized components, was investigated 
in paper I. Here, Figure 3 shows the location of both X (black 
symbols) and O (blue symbols) mode sources, as observed 
between 0812 and 0912 UT below f ce + 30 kHz, where both 
components were observed together. 

Either observed from Cassini (Figure 3a), or projected 
down to the planetary atmosphere (Figure 3b), distant X 
and O emissions originate from the same range of magnetic 
field lines, along the visible part of a spiral auroral oval, 
starting from very high latitudes (~ —80°) close to mid- 
night. Both components are mainly observed on the dawn 
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f-fce(kHz) f(kHz) 

Figure 4. (a) SKR beaming pattern of local sources (/ close to / ce ) over the interval [0812,0912] UT, 
expressed in beaming angle 0' — (k,-B) between the incoming wave direction and the magnetic field 
vector at the source, as a function of / — f ce . Each time- frequency k- vector direction is derived from 
GP analysis, with the data selection described in section 2, whereas the magnetic field vector at the 
source is either given by the SPV model and a simple current sheet (black symbols), or by in situ MAG 
measurements (green symbols). The median value computed over all measurements (solid line), and the 
standard deviation (dashed lines) are plotted for each distribution. SPV-derived 8' contains two distinct 
distributions [Cecconi el al, 2009]: reliable beaming angles are given by wave directions intercepting the 
surface / = f ce (crosses), while approximate beaming angles are estimated from wave that do not cross 
this surface (diamonds). Both distributions match within error bars with 8' ~ 70° below / CE + 1 kHz. 
Above this limit, MAG-derived 8' give less and less valid approximations of the beaming angle at the 
source and rather indicate the local wave propagation angle, (b) SKR beaming pattern of distant sources 
over the extended interval [0700,1100] UT and frequencies between 7 and 1000 kHz. Only SPV-derived 
9' are displayed here, separating LH (black symbols) and RH (blue symbols) polarizations, together with 
their median value (solid line) and standard deviation (dashed lines). In spite of scattered measurements 
(A9' ~ ±15°), the median beaming angle of the X mode (solid gray) displays a clear trend, with a 
plateau at 0'(f) ~75° below 300 kHz, followed by a decrease with frequency until #'(1000 kHz)~ 50°. 



side of the oval, between 0000 and 0800 LT. Nonetheless, 
a few X and O events also emanate from a dusk region at 
lower latitudes (~ —70°), possibly materializing the end of 
the spiral structure. Local sources (/ < f ce + 1 kHz, in red), 
identified as X mode in section 4, map the footprint of the 
field lines crossed by the spacecraft on Figure 3b. 

6. Beaming pattern 

The beaming pattern of planetary radio emissions is an- 
other crucial constraint to the generation mechanism at the 
source (when observed in situ), as well as to the wave propa- 
gation in the surrounding plasma (when observed remotely) . 
In this section, we compare the SKR beaming pattern de- 
termined in situ and remotely. 

Because investigated radio sources lie in the southern 
magnetic hemisphere, we will hereafter rather use 9' — 
n — 9 — (k,-B) to transpose results in the usual 0° to 90° 
range, and enable comparisons with previous work. 

6.1. Local sources 

We showed in section 4 that local sources are dominantly 
X mode emission. In Figure 4a, the corresponding distribu- 
tion of beaming angles 8' is derived over the interval of the 
source region by using both the SPV model (black) and in 
situ MAG observations (green), the latter being only valid 
locally. Results are plotted as a function of / — / ce , corre- 
sponding to gradually increasing distances d to the source. 

SPV-derived results consist of two distinct distributions: 
reliable beaming angles are given by wave directions inter- 
cepting the surface / = fee (crosses), while approximate 
beaming angles are estimated from wave directions that do 



not cross this surface (diamonds). Nonetheless, a reliable es- 
timate of 9' is given by the median computed over all beam- 
ing angle measurements (solid gray), as detailed in appendix 
A. The difference between the black and green distributions 
is related to the actual angle between the modeled field vec- 
tor at the source and the one measured in situ (the wave 
direction remaining unchanged). 

Below / ce + 1 kHz, both techniques tend to comparable 
results with 9' = 70° ± 15° for the black distribution and 
8' — 62° ± 6° for the green one (but with more and more 
dispersed measurements above 70° with decreasing / — f ce ). 
However, below / ce , the number of available measurements 
is not sufficient to give a reliable statistical estimate of 8' , 
that displays highly dispersed values, between 30° and 80°. 

To interpret this dispersion, we remind that both deter- 
minations of 8' rely on the validity of directions given by 
the 3-antenna GP inversion. Indeed, the correspondence be- 
tween black crosses and green triangles within a few degrees 
(that validates the use of the SPV model in the radio lo- 
calization technique) indicates that the angular uncertainty 
on the modeled field is negligible compared to the scatter- 
ing of k-vectors. This scattering is likely to be the result 
of two main factors. First, observations close to the source 
region can be affected by mixed emissions from several sep- 
arate intense sources. Second, as reminded in section 2, 
the accuracy of wave direction decreases with V, while Fig- 
ure lb specifically shows that V becomes low close to f ce . 
These reasons can consequently limit the accuracy of the 
determination of 8' for the closest sources, and we will keep 
8' = 70° ± 15° as a reasonable estimate of the beaming angle 
observed for local sources, quasi-perpendicular to the local 
magnetic field vector. 
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Above fce + l kHz, MAG-derived #'(/) shifts toward lower 
values, as the local magnetic field becomes a less and less 
valid approximation of the magnetic field at the source. In- 
stead, it simply indicates the local wave propagation angle. 
This frequency limit, already used in paper I and section 5, 
appears as a relevant upper limit for local sources. What- 
ever the frequency, the SPV-derived 9'(f) remains constant 
around ~ 75° ± 15°. 

In parallel with this method, based on the knowledge of 
the wave direction, an alternate technique based on polar- 
ization measurements only, and detailed in section 7, was 
used to measure reliable 9' for frequencies around f ce . We 
obtained 6' > 78° for local sources, that confirms and im- 
proves above results. 

6.2. Distant sources 

The SPV-derived beaming pattern of distant sources is in- 
vestigated in Figure 4b, along the extended time-frequency 
interval [0700,1100] UT and [7,1000] kHz, where emission is 
highly circularly polarized (Figure lb). 

LH (black symbols, V > 0) and RH (blue symbols, V < 0) 
polarized emissions were previously identified as X and O 
modes on the basis of the comparison of their intensity be- 
tween 15 and 30 kHz. At these frequencies, X and O mode 
components display similar 9' and comparable median val- 
ues (solid lines). However, close to 10 kHz, some RH events 
display intensities comparable to LH ones. As their esti- 
mated beaming angle is close to 90° (blue diamonds), they 
can be attributed to X mode, whose degree of circular po- 
larization changes sign at 9' = 90°. At frequencies higher 
than 40 kHz, where X mode prevails, the sign of V does not 
reverse at 90°, which suggests that beaming angles higher 
90° are not physical here and simply illustrate the statistical 
noise. 

We now focus on the X mode beaming angles. Extending 
results of Figure 4a, the median (gray solid) remains approx- 
imately constant at 9'(f) = 75° ±15° up to 300 kHz and then 
decreases with frequency toward #'(1000 kHz) = 50° ± 25°. 
In appendix A, we quantify with simulations the observa- 
tional bias induced by the employed radio localization tech- 
nique. In particular, we show that the median displayed 
here (computed over all LH measurements) gives a reliable 
estimate of real beaming angles, while the large scattering, 
increasing with frequency, can be attributed to the instru- 
mental noise on k directions and the observational geometry. 

In addition to the above direct determination, we also 
used an indirect technique to infer 9' for distant sources, that 
consists of modeling the narrowband intense part of the SKR 
low frequency envelope of Figure la, based on the PRES sim- 
ulation code (Planetary Radio Emissions Simulator) [Hess 
et al, 2008; Lamy et al., 2008b]. Here, we postulated SKR 
sources at frequencies covering the entire observed frequency 
range, along field lines located between 0030 and 0130 LT 
with invariant latitudes of —80°. As PRES assumes straight 
line propagation, 9' was then used a (sensitive) free param- 
eter to fit the low frequency envelope of the SKR spectrum 
between 10 and 60 kHz during the whole interval, as viewed 
mostly out of the sources. Best fit of the emission was ob- 
tained for 9' ~ 75°, with an angular width of the cone of a 
few degrees. This confirms SPV-derived 9' between 10 and 
60 kHz. 

These results can be compared to previous direct or indi- 
rect measurements of SKR beaming angles. Below 500 kHz, 
the median behavior of Figure 4b is marginally consistent 
with the 50° ±10° range derived for southern SKR sources by 
Cecconi et al. [2009] from 12 h of close mid-latitudes observa- 
tions. This could result from 9' variable with time, position 
and/or azimuth with respect to the magnetic field, as well as 
affected by refraction along the ray path, expected to play 



a non negligible role at low frequencies. A statistical study 
is requested to answer these questions. Interestingly, Fig- 
ure 4b is in excellent agreement with the oblique frequency- 
decreasing beaming angle computed by Lamy et al. [2008b] 
to model features in SKR dynamic spectra observed from 
equatorial latitudes at various distances. 

6.3. Intrinsic and apparent beaming pattern 

At Earth, AKR has been observed to be amplified quasi- 
perpendicularly to the local magnetic field within its source 
region [Hilgers et al., 1992], from unstable trapped [Louarn 
et al, 1990], or shell- type [Ergun et al., 2000] electron dis- 
tributions. Then, refraction along the path, at the bound- 
aries of auroral cavities or through the terrestrial plasmas- 
phere, results in a complex oblique cone of emission, non- 
axisymmetric and partially filled [Louarn and Le Queau, 
1996b, a; Xiao et al, 2007; Mutel et al., 2008]. 

Here, local sources below / CE ±1 kHz display large 9' , mea- 
sured between 30° and 80° from direction measurements, or 
above 78° from polarization measurements, consistent with 
perpendicular emission within error bars. This result is sup- 
ported by conclusions drawn in companion studies: in paper 
I, we inferred perpendicular emission from the comparison of 
SKR cutoff frequencies and CMI resonance frequencies. In 
parallel, Mutel et al. [2010] computed CMI growth rates us- 
ing ring-type electron distributions fitted to shell-like struc- 
tures observed within the source region [Schippers et al., 
2010], that account for the observed SKR intensity by being 
maximal for perpendicular amplification. In addition, we 
can notice that 9' > 70° together with a cone width of 5°, 
as previously estimated from far observations [Lamy et al., 
2008b], yields a solid angle of ~ 0.5 sr, anticipated in section 
4.4. 

Above f ce ± 1 kHz, or equivalently for distances beyond 
~ 0.5 Rs, beaming angles remain constant at 9' — 75° ±15°, 
considered as a reliable estimate, up to 300 kHz, then de- 
creasing with frequency. Lamy et al. [2008b] modeled such 
a frequency-dependent beaming angle from a loss cone elec- 
tron distribution. As CAPS measurements rather suggest 
shell or ring-type distributions, that favor perpendicular 
emission, alternate possibilities must be explored. Emit- 
ted radio waves are likely to undergo refraction within a few 
tenths of kronian radii until they eventually reach 9' ~ 75° 
in regions where N ~ 1. Indeed, in Figure 3, the locus of lo- 
cal sources precisely matches the footprints of magnetic field 
lines that were crossed by Cassini while assuming straight 
line propagation. This confirms that refraction must occur 
very close to the amplification region. We can make sim- 
ple estimates of refractive conditions along the observed ray 
paths: assuming 9' = 75° with N — 1 for distant sources, 
the Snell-Descartes law requires N ~ 0.966 at the source for 
a strictly perpendicular emission. 

Finally, the striking similarity of X and O mode beaming 
angles is intriguing. If both modes emanate from the same 
regions, as suggested by Figure 3, at frequencies close to f ce , 
itself much larger than f pe , refraction is expected to be more 
significant for the X mode close to fx (N x < 1), than for 
the O mode far above fo (N () ~ 1). 

A dedicated ray tracing study is required to quantify the 
SKR ray path over all accessible directions for accessible 
magneto-ionic modes. Interestingly, we note that the jovian 
decametric emission (DAM) also displays beaming angles 
about 70° [Queinnec and Zarka, 1998]. 

7. State of polarization 

7.1. SKR polarization ellipse along the propagation 

Previous studies showed that SKR is quasi-fully polar- 
ized, with a degree of total polarization (V 2 + Q 2 + U 2 ) 1 ^ 2 ~ 
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Figure 5. Polarization state of radio waves observed between 0712 and 0912 UT (source region extended 
by one hour to include significant O mode emission) and f ce —2 kHz < / < f ce +30 kHz. (a) Circular and 

(b) linear degree of polarization V and L = (Q 2 + U 2 ) 1//2 are displayed as a function of / — f ce , whereas 

(c) / - fee and (d) axial ratio of the polarization ellipse T = [(Q 2 +U 2 +V 2 ) 1/2 - L]/V [Kraus, 1966] 
are displayed as a function of the distance to the source (the distance along which the wave direction 
intercepts the surface / = f CB ). Black and blue symbols refer to V > and V < respectively identified 
as X and O modes in section 4.1. As these modes are expected to be orthogonal (T X T D = — 1), Figure 
5d superimposes T x and 1/T D , so that polarization ellipses can be directly compared. RPWS data 
have again been selected with the selection described in section 2, with an additional criterion retaining 
only measurements with a normalized degree of total polarization (Q 2 +U 2 +V 2 ) 1//2 between 0.85 and 1.1. 
Moreover, considering that local sources correspond to / < f ce + 1 kHz, as illustrated in Figures 3 and 
4, associated distances have been systematically set to 0, even when not crossing the surface / = f ce , X 
and O modes follow a similar trend in all panels and d regularly increases with frequency, as expected. 
SKR waves are observed to be quasi-purely circularly polarized above fee + 10 kHz (d > 1.5 Rs), and 
elliptical below, toward V ~ 0.2 - 0.4, L ~ 0.8 - 0.9 and T x ~ 0.1 - 0.3 at / = /„. 



1, and that the observed polarization evolves with the loca- 
tion of the observer [Fischer et at, 2009, and references 
therein] . Whereas SKR polarization is quasi-purely circular 
when detected from low latitudes (< 30°), it becomes mainly 
and significantly elliptical when observed from higher lati- 
tudes (between 30° and 60°). 

Here, Cassini's trajectory swept very high southern lat- 
itudes, from —40° to —75°, along which two important 
features are brought to light by Figure lb. First, distant 
sources above f ce display high levels of circular polarization 
with | V| ~ 1 (in regions of the spectrogram where GP results 
are not affected by the geometric configuration between an- 
tennas and real wave directions). This strongly differs from 
the trend noticed by Fischer et al. [2009] in the range —40° 
to —60°. Second, as noticed in section 4.1, circular polar- 
ization vanishes for local sources, close to f ce , whereas we 
checked that radio waves remain quasi-fully polarized. This 



suggests that the wave polarization changes along the ray 
path with the distance of propagation. 

Figure 5 investigates the SKR full state of polarization in 
detail. As mentioned earlier, 3-antenna measurements bring 
a direct tridimensional determination of the polarization el- 
lipse. On the left-hand side, the degrees of circular (V, Fig- 
ure 5a) and linear (L, Figure 5b) polarization are plotted as 
a function of / ' — f ce . On the right-hand side, f — f ce (Figure 
5c) and the axial ratio of the polarization ellipse (T, Figure 
5d), are plotted as a function of the distance to the source 
d. Figure 5c illustrates that the quantity / — f ce is a proxy 
for d. The dataset employed here retains only quasi-fully 
polarized events with (V 2 + Q 2 + U 2 ) 1/2 > 0.85, observed 
between 0712 and 0912 UT (to include significant X and O 
mode contributions) , and frequencies below f ce + 30 kHz (d 
between and 3 Rs)- 

SKR polarization continuously evolves with / — f ce . 
Above 10 kHz, the observed radio sources are located be- 
yond 1.5 Rs and are highly circularly polarized, with |V| ~ 
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0.8 - 1.1, L ~ 0.0 - 0.3 and T ~ 0.7 — 1.0. A few events, 
seen at all frequencies with low V , high L and low T, cor- 
respond to unreliable GP results (located in time-frequency 
regions of weak signal on Figure lb). Below 10 kHz to- 
ward decreasing / — fee, a continuous trend shows more and 
more elliptically polarized emission, reaching V ~ 0.2 — 0.4, 
L ~ 0.8-0.9 and AR ~ 0.1-0.3 at the limit / = f ce . Stokes 
parameters are very similar for X and O modes, whatever 
the frequency. This is not surprising since characteristic po- 
larizations are expected to be orthogonal (T X T G = —1, see 
section 3), which means that polarization ellipses of X and 
O mode only differ by a rotation of n/2, the former having 
its main axis along the magnetic field direction. 

In summary, SKR polarization is naturally strongly ellip- 
tical at the source, before being circularized along its prop- 
agation path. These two aspects are investigated below. 

7.2. Intrinsic polarization 

In contrast to the common picture of planetary auroral 
radio emission being circularly polarized, SKR naturally dis- 
plays elliptical polarization at the source. How does this 
compare to predictions of the magneto-ionic theory? 
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Figure 6. Identical to Figure 5d, except that axial ra- 
tios predicted by equation 3, using X , Y and local wave 
propagation angles 8' , have been overplotted with green 
triangles. The overall correspondence between both dis- 
tributions is clear. 




Figure 7. Coupling factor Q [Sawyer et al, 1991], com- 
puted from equation 6 and data of Figure 5d, as a func- 
tion of d. Q remains small whatever the distance of prop- 
agation. 



Using conditions prevailing within a CMI source (X <C 1, 
Y ~ 1 and weakly relativistic electrons), Melrose and Dulk 
[1991] showed that the polarization of a cyclotron emission 
and the one of the X mode coincide and reduce to T ~ cos 8' 
to first order. Here, we can use the exact expression given by 
equation 3 with local measurements of T x , X and Y to di- 
rectly compute the beaming angle at the source. For closest 
sources in Figure 5d, T x < 0.2 yield 8' > 78°, supporting 
quasi-perpendicular emission suggested in section 6. 

The relation T~ cos 8' is not valid for distant sources, 
for which Y shifts from unity. As an illustration, while ax- 
ial ratios increase with d (Figure 5d), 8' remains constant 
(Figure 4). Nonetheless, the knowledge of X, Y as well as 
the local wave propagation angle (green 8' in Figure 4) for 
each time-frequency measurement allows us to directly use 
equation 3 to predict the local value of T x (f) and 1/T G (/) 
for all waves (including local and distant sources), whatever 
their ray path. 

Observed and predicted axial ratios are superimposed in 
Figure 6. Though some waves propagate through distances 
of a few kronian radii, the overall correspondence is pretty 
good. Precisely, T x (f) and 1/T Q (/) evolve similarly be- 
tween 0.5 and 2 Rs- Below 0.5 Rs, predicted T x are slightly 
higher than observed ones, which is the consequence of un- 
derestimated MAG-derived 0'(f) for local sources (~ 70°). 
Above 2.5 R cf: , the discrepancy between observed and pre- 
dicted T x corresponds to noisy results observed close to 
40 kHz in Figure lb. In summary, the magneto-ionic the- 
ory satisfactorily accounts for the observed SKR elliptical 
polarization. 

7.3. Polarization transfer 

The above result can be analyzed in the frame of the 
mode coupling theory [Booker, 1936; Budden, 1952; Cohen, 
1960; Melrose, 1980; Daigne and Ortega- Molina, 1984], used 
to investigate the elliptical polarization of jovian decamet- 
ric emission (DAM) [Goertz, 1974; Lecacheux et al, 1991; 
Melrose and Dulk, 1991; Shaposhnikov et al, 1997] and pro- 
posed to analyze first results on SKR elliptical polarization 
[Fischer et al, 2009]. 

When propagating through an inhomogeneous plasma, 
the polarization transfer of a radio wave can be described 
by two limiting regimes. If the medium is changing slowly 
enough (X, Y and 8' do not change significantly over a 
wavelength), it can be considered as locally homogeneous 
and characteristic magneto-ionic modes are "weakly cou- 
pled" . This means that natural modes are preserved and 
propagate independently of each other. The magneto-ionic 
theory prevails and their polarization and index of refrac- 
tion, described by equations 1-3, slowly evolve along the ray 
path. For instance, their sense of circular polarization can 
reverse when k.B changes sign (8 — tt/2 in equation 3), 
and Faraday rotation can occur. In contrast, if the medium 
changes rapidly, characteristic modes become "coupled" and 
the magneto-ionic theory does not apply any more. Char- 
acteristic modes do not propagate independently and can 
transfer energy to each other, so that the polarization of a 
given wave is retained along the ray path. The radiation 
propagates as in vacuum. 

The transition between regions of weak and strong cou- 
pling has been interpreted as a limiting polarization zone 
(LPZ), beyond which the wave polarization is fixed [Budden, 
1952] . It has been characterized by Cohen [1960] in the high 
frequency approximation (X, Y <C 1) in terms of a transi- 
tional frequency ft whose expression, established in the cases 
of quasi-longitudinal (QL) and quasi-transverse (QT) prop- 
agation, is proportional to the plasma density and the mag- 
netic field strength (see precise expressions in equations 53 
and 58). Cohen [1960] showed that weak coupling occur for 
/ <C ft or f 4 <C f 4 respectively for QL or QT propagation, 
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and strong coupling when the opposite condition is fulfilled. 
For the jovian magnetosphere, Lecacheux [1988] reduced ft 
to a function of the electron density n, and showed that 
strong coupling at 30 MHz requires n < 5 cm . A similar 
result was obtained by Melrose and Dulk [1991] for a region 
close to the source, using Y ~ 1/2 and a typical length of 
inhomogeneity of one jovian radius, and then adapted to 
Saturn by Fischer et al. [2009] to infer n < 0.01 cm~ 3 for 
SKR frequencies below 1 MHz. Transposed to the equiv- 
alent frequency of 20 kHz in the present study (Y ~ 1/2 
and fee ~ 10 kHz), the region of strong coupling shall be 
reached for n < 0.001 cm , a much lower value than the 
0.01 cm~ 3 measured locally (see paper I). This qualitatively 
shows that the LPZ at 20 kHz has not been reached yet and 
that associated radio waves propagate in a region of weak 
coupling. In addition, it is unlikely to find such low elec- 
tron densities within the kronian magnetosphere, so that 
the LPZ may not be reached until waves eventually go out 
of the magnetosphere. 

To estimate precisely the amount of coupling along the 
ray path, it is more convenient to use the coupling factor 
Q introduced by Budden [1952], and defined as the ratio of 
the rate at which the polarization of one characteristic mode 
changes along the ray path to the rate at which character- 
istic modes get out of phase. For instance, Sawyer et al. 
[1991] explicited Q as: 



1 dT/dd{T 2 - I)-' 1 
Q ~2 2vf\AN x J/c 



(6) 



where AN X = \N Q — N x \. Regions of weak and strong 
coupling then correspond to Q <JC 1 and Q> 1, respectively. 
Here, dT/dd can be simply estimated from Figure 6a, where 
T x and 1/T vary approximately linearly by an amount of 
0.7 over a distance of 2 Rg, while AN xo is computed from 
equation 1 from X, Y and 9' . 

Figure 7 shows Q(d) computed for X and O mode emis- 
sions displayed in Figures 5d and 6. This result was checked 
by using alternate coupling factors derived in other polariza- 
tion studies [Daigne and Ortega-Molina, 1984; Melrose and 
Dulk, 1991; Shaposhnikov et al, 1997], that all give similar 
or lower values of Q. Low coupling factors are thus observed 
at all frequencies, with d < 3 Rs, and consequently charac- 
terize a region where characteristic modes are not or weakly 
coupled. SKR is mainly emitted and remains in the X mode, 
with characteristics described by the magneto-ionic theory. 
An important consequence is that Faraday rotation is un- 
likely to occur along the investigated ray paths. 

7.4. Consequences and expectations for other planetary 
auroral radio emissions 

Among planetary auroral radio emissions, elliptical po- 
larization was unambiguously observed only for jovian DAM 
and kronian SKR, whereas AKR at Earth or UKR at Uranus 
where characterized as being mostly circular, as reminded by 
Fischer et al. [2009]. 

At Jupiter, early DAM studies proposed that the ob- 
served elliptical polarization could be fixed by strong cou- 
pling at the LPZ in the close neighborhood of the source 
[Lecacheux et al., 1991; Melrose and Dulk, 1991], in which 
case the wave polarization would carry information on the 
source region and the generation mechanism. Alternately, 
Shaposhnikov et al. [1997] proposed that elliptical polariza- 
tion might be formed by moderate linear mode coupling in 
transitional regions (TR) outside of the source region, de- 
fined by A 2 ~ 1 at the transition between QT and QL prop- 
agation, and carries little information on the source region. 
Here, X = 0.01, 9' = 75° and Y varying from 1 to 1/2 corre- 
sponds to A 2 between 3.3 and 0.8. Thus, when propagating 
out of their source region, SKR waves precisely lies within 
the transitional region where weak mode coupling occurs. 

At Earth, quantitative measurements of AKR polariza- 
tion with several spacecraft have established that the emis- 
sion is strongly circularly polarized, whatever the location of 



the observer [Lefeuvre et al, 1998; Panchenko et al., 2008; 
Larny et al., 2010b]. At the light of SKR polarization proper- 
ties, the non detection of elliptical AKR is intriguing. From 
the expression of the transitional frequency transposed to 
Earth, weak coupling is expected for n > 0.1 cm~ 3 at the 
AKR frequency peak of 200 kHz, a condition which is ful- 
filled both within (n ~ 1 cm~ 3 ) and outside (n ~ 5 to 
10 cm ) the AKR source region [Roux et al, 1993; Ergun 
et al, 1998]. We consequently expect AKR polarization to 
vary as the one of the characteristic mode on which it is 
emitted. Owing to similar conditions of generation for SKR 
and AKR (large 9, Y ~ 1 and X <C 1), the AKR polar- 
ization at the source should be strongly elliptical at least 
within auroral cavities. 

At Uranus, a polarization model based on the magneto- 
ionic theory was developed to derive the limiting polariza- 
tion of Uranian Kilometric Radiation (UKR) [Sawyer et al., 
1991], and successfully retrieved axial ratios derived from a 
few 2-antenna measurements of Voyager. This strengthens 
the picture of planetary radio emissions polarized accord- 
ingly to the conditions encountered along the ray path. 

Coming back to Saturn, considering that SKR emission 
is beamed in a thin cone of emission aligned with the lo- 
cal magnetic vector, the plasma parameters X , Y and the 
wave propagation angle 9' will change differently with the 
azimuthal angle around B (in spherical coordinates) , as well 
as the source locus. For instance, a radio wave propagating 
toward the equatorial dense plasma will encounter increas- 
ing densities, preventing strong coupling and circularizing 
the wave polarization, in agreement with near-equatorial 
observations below 30° latitude. As the index of refraction 
and the wave polarization are expected to vary significantly 
with azimuth, invariant latitude and local time, a detailed 
modeling of the SKR polarization transfer in appropriate 
plasma conditions, beyond the scope of this paper, is re- 
quired to determine why elliptical polarization is well de- 
tected over the whole SKR spectrum for mid-latitudes ob- 
servations and which information on the source is carried by 
observed waves. 



8. Summary 

In this study, we measured the properties of SKR along 
its propagation path from its source region to the observer 
thanks to simultaneous radio (electric and magnetic wave 
components), magnetic field and electron measurements ac- 
quired by the Cassini spacecraft. 

SKR is mainly emitted in the fast X mode, at frequen- 
cies above the X mode cutoff derived for the observed hot 
auroral plasma, dominated by 6 to 9 keV electrons. From in 
situ observations, we estimated the mean electron-to-wave 
energy conversion efficiency about 1% (2% peak), in agree- 
ment with previous macroscopic estimations. SKR is also 
observed in the O mode, typically 2 orders of magnitude 
fainter, whereas no kilometric Z mode was detected. 

X and O SKR components are found to emanate from 
the same regions, and display similar beaming pattern and 
polarization at comparable frequencies. The SKR beaming 
pattern has been estimated by various techniques. The ra- 
dio localization technique allows us to measure a rough cone 
of emission with large aperture angle 9' — 70° ± 15° below 
fee + 1 kHz, corresponding to the upper frequency limit for 
local sources. But the most accurate beaming angles are 
given by direct polarization measurements yielding 9' > 78°, 
which supports a quasi-perpendicular emission, consistently 
with results inferred in companion studies. For distant 
sources above f ce + 1 kHz, 9'(f) = 75° ± 15° up to 300 kHz 
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Figure 8. (a) Southern polar view of a simple spiral auroral oval, comparable to the one of Figure 3b. 
(b) Equatorial view from dawn of corresponding magnetic field lines (gray) , with simulated radio sources 
distributed between 6 and 1000 kHz (black crosses). The location of the observer is set at 08:12 (green 
triangle) along the Cassini's trajectory (green dashed). Blue crosses display radio sources visible from 
the spacecraft for a beaming angle 9' (here equal to 75° ±5°). 



and decreases above to reach 6' (1000 kHz) = 50° ±25°. The 
rapid transition from quasi-perpendicular to oblique wave 
propagation angles could result from refraction along the 
ray path close to the source region. These results signifi- 
cantly differ from previous measurements of SKR beaming 
angles, likely to result from beaming angles varying with the 
location of the observer. 

SKR polarization is shown to evolve with propagation, 
from strongly elliptical at the source to highly circular af- 
ter propagation across 2 Rs of high-latitude auroral plasma. 
Thanks to the knowledge of the local plasma parameters and 
wave propagation angle, we demonstrated that the polariza- 
tion predicted by the magneto-ionic theory satisfactorily de- 
scribes the observations. Conditions of weak mode coupling 
apply to the present event from the sources to a distance 
larger than 3 Rs, so that the SKR polarization corresponds 
to the one of characteristic modes. In addition, similarly to 
SKR, we predict that AKR polarization measured (at least) 
inside its source region shall be strongly elliptical. 

Considering the numerous striking similarities between 
AKR and SKR, measured remotely and locally, the CMI 
process, responsible for AKR generation is reinforced as a 
universal mechanism responsible of planetary auroral radio 
emissions. The differences noticed between AKR and SKR 
source regions (frequency of emission, presence of auroral 
cavities, acceleration processes, observed polarization) are 
characteristic of the CMI operation in each magnetosphere. 

Appendix A: Simulation of radio source 
beaming pattern and observational bias 

The determination of beaming angles from the radio lo- 
calization technique used in this article has been described 
and discussed by Cecconi et al. [2009]. For radio waves ob- 
served at a frequency /, reliable 9' are obtained from direc- 
tions that intercept the surface / = f ce (crosses in Figure 
4), whereas only approximate 9' (generally close to 90°) can 
be derived from directions that do not intercept this surface 
(diamonds in Figure 4). These two distributions are gen- 
erally distinct, as a result of the geometrical configuration 
between visible sources and the observer. 

Here, we want to determine how Figure 4b is affected by 
this observational bias, and check if the median computed 



over all beaming angle measurements yields a reliable esti- 
mation of the real beaming angle. 

To do so, we build a simple spiral auroral oval (Figure 
8a), comparable to the one of Figure 3b. Along magnetic 
field lines mapping to this oval, we simulate southern radio 
sources distributed between 6 and 1000 kHz (Figure 8b). 
Then, we set the spacecraft location at 08:12 (green trian- 
gle), corresponding to the entrance into the SKR source re- 
gion, and identify which radio sources are visible for a fixed 
beaming angle 9' (blue crosses, here using 9' = 75° ± 5°) 
with a straight line propagation (N = 1). 

Afterwards, we model errors on visible sources by adding 
a gaussian noise distribution, made of 100 points and an 
angular width at half maximum equal to 2° [Cecconi and 
Zarka, 2005], to each direction of arrival (Figure 9a, c). This 
blurred dataset is processed by the radio localization tech- 
nique, and finally provides simulated beaming angles that 
directly compare to 9' (Figure 9b, d). 

Figures 9a, b and 9c,d display results obtained for 9' = 
75° ± 5° and 9' — 90° ± 5°, respectively. As expected, 
visible radio sources lie on the dawnside, with a blurred dis- 
tribution comparable to the one of Figure 3a, whereas slight 
differences of visibility are related to the value of 9' . Our 
purpose here is not to strictly reproduce observations but to 
simply approach a realistic situation. 

The distributions of simulated beaming angles allow us 
to quantify the geometrical bias for observing conditions of 
Figure 4. First, the distribution of reliable beaming angles 
(crosses) naturally enlarges with frequency as a result of 
the distance to the source, with a median (gray dotted) de- 
creasing with frequency. The latter artificially shifts from 9' , 
which is only reached only at the lowest frequencies. Second, 
the distribution of approximate beaming angles (diamonds) 
is close to 90° and slightly increases with frequency. Both 
of these features are well observed on Figure 4b. 

Then, we note that the median computed over all beam- 
ing angles (solid gray) remains approximately constant and 
equal to 9' whatever the frequency and the value of 9' , which 
means that it provides a reliable estimation of real beaming 
angles. As a result, we conclude that the median displayed 
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Figure 9. (a) Blurred distribution of radio sources as viewed from the spacecraft with 8' = 75° ± 5°, 
comparable to Figure 3a. (b) Corresponding distribution of beaming angles, once processed by the radio 
localization technique. Similarly to Figure 4, crosses (diamonds) refer to directions of arrival that cross 
(do not cross) the surface / = f ce . The dashed gray line corresponds to the median computed over reli- 
able beaming angles (crosses) only, while the solid gray line gives the median computed over all beaming 
angles (crosses and diamonds). The latter is approximately equal to 8' , whatever the frequency. (c,d) 
identical to (a,b) with 0' = 90° ± 5°. 



in Figure 4b (approximately constant around ~ 75° up to 
~300 kHz, decreasing at higher frequencies) cannot be ex- 
plained by a geometrical bias and shows a real physical ef- 
fect, either due to the driving of the CMI by a loss-cone 
electron distribution [Hess et al, 2008; Lamy et al, 2008b] 
or to frequency-dependent refraction. 
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